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M84 - A Warp Caused By Jet Induced Pressure Gradients? 1 

A. C. Quillen 2 & Gary A. Bower 3 

ABSTRACT 

| In radio galaxies such as M84 dust features tend to be nearly perpendicular to 

radio jets yet are not aligned with the galaxy isophotes. The timescale for precession in 
the galaxy is short (~ 10 years at 100 pc) suggesting that an alternative mechanism 
causes the gas disk to be misaligned with the galaxy. In M84 we estimate the pressure 
on the disk required to overcome the torque from the galaxy and find that it is 
small compared to the thermal pressure in the hot ambient ISM estimated from the 
X-ray emission. We therefore propose that pressure gradients in a jet associated hot 
interstellar medium exert a torque on the gas disk in M84 causing it to be misaligned 
' with the galaxy isophotal major axis. We propose that AGN associated outflows or 

associated hot low density media in nearby radio galaxies could strongly affect the 

■ orientation of their gas disks on 100 pc scales. This mechanism could explain the 
connection between gas disk angular momentum and jet axes in nearby radio galaxies. 

, By integrating the light of the galaxy through a warped gas and dust disk we find 

\ that the geometry of gas disk in M84 is likely to differ that predicted from a simple 

precession model. We find that the morphology of the gas disk in M84 is consistent 
with a warped geometry where precession is caused by a combination of a galactic 

■ torque and a larger torque due to pressure gradients in the ambient X-ray emitting 
gas. Precession occurs at an axis between the jet and galaxy major axis, but nearer 

■ to the jet axis implying that the pressure torque is 2-4 times larger than the galactic 
torque. We estimate that precession has occurred about this particular axis for about 
10 7 years. A better model to the morphology of the disk is achieved when precession 
takes place about an elliptical rather than circular path. This suggests that the isobars 
in the hot medium are strongly dependent on angle from the jet axis. 
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1. Introduction 



Recent imaging from the Hubble Space Telescope (HST) and spectroscopic studies of active 
elliptical galaxies have established that there are gas disks in the central few hundred pc of these 
galaxies which could be feeding the massive ~ 10 9 M Q black holes at their centers (e.g M84; Bower 



et al. 1997a, Bower et al. 1997b, M87; Harms et al. 1994 and NGC 4261; Ferrarese, Ford, fc Jaffe 



1996). These gas disks are often observed to be nearly perpendicular to the jets, establishing a link 



between the angular momentum of the disk and the direction of the jet. Surveys of the observed 
angular difference, between radio and disk axes in radio ellipticals and find a statistically 
significant peak in the distribution at \l/ = 90° ( Kotanyi fc Ekers 1979| , Mollenhoff, Hummel & 



Bender 1992). This has established that radio jets are generally perpendicular to gas disks on 



large (kpc) scales in nearby radio galaxies, even though no correlation between the radio axis and 
the galaxy isophotal major axis exists at low redshift (e.g., Sansom et al. 1987| ) . The strong peak 
at = 90° has been confirmed at smaller (100 pc) scales by van Dokkum fe Franx (19951 
HST/WFPC images. 



using 



Strangely, the jet direction is not necessarily expected to be correlated with the angular 
momentum of the gas feeding the black hole (and so the orientation of the disk) because of 
Lense-Thirring precession, otherwise known as 'dragging of inertial frames' flRees 19"78 ). A 



spinning black hole causes a gas disk near the gravitational radius, 



GM/c 2 



10 



ii 



cm, 



to precess, and so the disk or torus near the black hole is expected to fill a volume that is 
axisymmetric and aligned with the spin axis of the black hole itself ( [Rees 19*84 , Bardeen & 



Petterson 1975). This inner torus is the proposed site for jet collimation and acceleration (Rees 
et al. 1982] ) so we expect that the jet should be aligned with the spin axis of the black hole but 
not necessarily the angular momentum of the disk well outside of r g . A massive spinning black 
hole can be described as an angular momentum reservoir that does vary in momentum but only 
very slowly with the influx of fuel. The rate of change of angular momentum is expected to be 
particularly slow (;> 10 9 years) in radio galaxies where the black holes are massive (~ 1O 9 M0) and 
the fueling rates are probably low (~ lO _4 M0/yr; [Rees 1978] ), though accretion of a second black 
hole could change the spin axis of the AGN on a very short timescale. If, however, the black hole 
does not spin significantly, then the jet axis would be determined by the angular momentum of 
the inner disk. This disk, because of its lower mass, could vary in orientation on shorter timescales 
than the black hole axis (for a spinning system). However, the observed alignment of jet and disk 
axes would then require that the orientation of this inner disk be coupled to that of the outer disk 
(at ^ lOOpc). Possible mechanisms causing alignment between the jets and a disk at 100 pc scales 
(well beyond r g ) have not been explored. 

To explore alignment mechanisms, we consider the geometry of the warped dust features 
and ionized emission in M84 (NGC 4374, 3C 272.1), a radio elliptical in the Virgo cluster. M84 
has a dusty disk seen in extinction in optical images (see Fig. 1) and in emission in Ha +[NII] 



that although nearly perpendicular to the jet, is misaligned with the galaxy isophotes ( Hansen 



Norgaard-Nielsen, fc Jorgensen 1985, van Dokkum fc Franx 1995, Jaffe et al. 1994, Bower et 
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al. 1997a, Baum et aL 1988] ). |van Dokkum fc Franx (1995) noted that the timescale for the disk to 
settle into the galaxy plane of symmetry is probably short; ~ 10 s years. These authors therefore 
suggested that the misalignment could be explained by a rapid inflow of gas into the central region 
of the galaxy (to radii smaller than r <^ 600 pc). However if this is the case, it is difficult to 
explain the S-shape of the Ha + [Nil] emission ( [Baum et al. 1988| ) and the detailed morphology 
of the dusty disk itself. This emission is elongated and lies within 20° of PA ~ 70° for r <^ 7" but 
at this radius the disk twists to PA ~ 115° on either side of the nucleus. It is difficult to imagine 
how a cooling flow could generate such a morphology. As an alternative to this possibility, we also 
consider an interaction between the gas disk and a hot ambient interstellar medium (ISM) as a 
possible mechanism for affecting the alignment of the disk. By integrating the light of the galaxy 
through a warped gas and dust disk we find that the geometry of gas disk in M84 is likely to differ 
that predicted from a simple precession model. We search for a simple physical model that can 
account for the observed disk morphology. 

Throughout this paper we adopt a distance to M84 of 17Mpc flMould et al. 1995| ; H = 75km 
s~ 1 Mpc^ 1 ). At this distance 1" corresponds to 82 pc. 



2. Misalignment of the Dust Features with the Galaxy 



That the dust features are misaligned with the galaxy isophotes in M84 for r <J 600 pc 
was emphasized by van Dokkum &: Franx (1995)| . The galaxy isophotes have a major axis PA 
= 129° ( |van Dokkum fc Franx 1~995| ; see Fig. 1) whereas the dust features are elongated along 
PA= 80 — 65° (for r < 7") depending upon the region used to estimate the angle ( Bower et 



al. 1997a , van Dokkum fc Franx 1995 ). This is a misalignment of 50 
major axis with respect to the galaxy major axis. 



65° of the dust features 



Misalignment between gas and stellar isophotes is not necessarily rare, (e.g. see van Dokkum 



fc Franx 1995 ). When there are kpc scale gas/galaxy misalignments usually the misalignment is 
the result of a merger (e.g. Tubbs 1980 and |Steiman-Cameron, Kormendy, fc Durisen 1992; ). A gas 
disk in a non-spherical but axisymmetric galaxy will precess about the galaxy axis of symmetry. 
At a radius of a few kpc the time for a gas disk to precess about the galaxy axis of symmetry 
is long, ~ 10 8 years, requiring a merger event to have occurred within this time if a merger is 
responsible for the gas/galaxy misalignment. Because the precession is faster in the center of the 
galaxy in the central regions the warp can be multiply folded and in the outer regions where the 
precession is slower the angle of the disk is more directly related to the orbital angular momentum 
of the merger event ( |Quillen, Graham, fc Frogel 1993] ). A multiply warped disk will form a band of 
absorption that is roughly perpendicular to the galaxy axis of symmetry. The morphology of the 
dust lanes in M84 resembles those of Centaurus A and NGC 4753 which are caused by multiply 
folded thin dusty surfaces ( |Quillen et al. 1993| , |Steiman-Cameron et al. 1992 ). However precession 
in an axisymmetric galaxy potential is an unlikely explanation for the warp in M84 since the band 
of absorption in the multiply folded region is not aligned with the galaxy isophotes. 
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van Dokkum & Franx (1995) emphasized that the timescale for the gas disks to settle into the 
plane of the galaxy is short in the central few hundred pc. The timescale for a gas ring to precess 
about the galaxy axis of symmetry, is shorter still. The angular precession frequency da/dt about 
the galaxy axis of symmetry for a gas ring of radius r and inclination 6 g with respect to the galaxy 
axis of symmetry is 

do. 

— = e$f2 cos 9 g (1) 
(e.g. Gunn 1979] ) where e$ is the ellipticity of the gravitational potential which we assume is 



axisymmetric and Q = v c /r for v c the velocity of a particle which has only a gravitational force on 
it in a circular orbit. Except in the case of a polar ring {0 g ~ 90°), the dependence on inclination 
is weak. In a time 

which we refer to as a precession time, rings that differ in radius by a factor of 2 should have an 
angular difference of ~ tt/2. 

To calculate the precession timescale we estimate the circular velocity, v c ~ l.lcr* (for a 
mildly anisotropic velocity distribution) where cr* = 303 ± 5km/s for r <^ 5" ( Davies fc Birkinshaw 



1988[) is observed stellar velocity dispersion. The ellipticity of the gravitational potential, e$, 



which is directly related to the ellipticity of the isophotes (see Fig. 1) e$ ~ e/3 (see |Binney & 



Tremaine 1987 , Figure 2-13). Here the isophote ellipticity, e = 1 — b/a, where b/a is the axis ratio 
of the isophotes. In the near-infrared NICMOS 2.2/xm and 1.6/um images the isophotes of M84 are 
symmetrical, do not vary in position angle, and show little extinction from dust. The ellipticity 
from these NICMOS images shows that an e > 0.1 persists to well within the core break radius 
(?-{, ~ 2") to r ^ 0.3" ( Quillen, Bower Stritzinger 1999) ) . This results in 



Since in M84 the dust is misaligned with the galaxy for r <J 600 pc, a mechanism operating at a 
timescale faster than the above timescale must be causing the disk to remain at an angle differing 
from the galaxy isophotes. 



2.1. Pressure Required to Overcome the Galaxy Torque 



Here we consider the possibility that the misalignment of the disk with the galaxy is caused 
locally. Radio galaxies put a substantial fraction of their luminosity into the kinetic energy of their 
jets (e.g. Rees et al. 1982|) . Emission in Ha + [Nil] commonly shows high velocity dispersions of 
the order of a few hundred km/s ( |Baum, Heckman, fc van Breugel 1992 , Tadhunter, Fosbury, & 
Quinn 1989, Axon et al. 1989]) . In radio galaxies there is strong evidence that the jets themselves 



impart significant motions in the surrounding ISM, implying that a significant fraction of the jet 
mechanical energy is dissipated locally in the surrounding ISM (e.g. De Young 1981; Begelman 
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1982). If the jets are responsible for motions in the surrounding ISM then there should be a 
differential in this medium, with the largest pressures nearest the jets, and the lowest pressures in 
the plane perpendicular to the jets. 

If the disk orientation is affected by the local medium then there must be a pressure on the 
gas disk greater than that of the gravitational potential which would be causing it to precess. The 
torque per unit mass on a ring of radius r is 

t = 6$ cos (Og) sin (9g)v^. (4) 

and so the pressure required to keep the disk from precessing is 

P T > 6$ cos (fig) sin (6 g )Ev* /r (5) 

where £ is the mass per unit area in the gas disk. This pressure is 



_ 2 / e$ \ / E \ / v c \ 2 /lOOpc 



P r -Z 2 * 10-dynes cm- ) (jj^) (jj^j ) « (»,) sin (0,) (6) 

and can be compared to that estimated from the hot medium based on the X-ray emission. From 



the X-ray emission, Thomas et al. (1986) find an electron density that is n e ~ 0.5(100pc/r)cm -3 
(where we have extrapolated from their profile which ranges from 0.5 - 20 kpc). This extrapolation 
is somewhat justified by comparison of pressures estimated from the [SII] 6717, 673lA lines in 
cooling flow galaxies (including M87 within r < 200 pc) which are not inconsistent with high 
central densities extrapolated from the lower resolution X-ray estimated pressures ( [Heckman et 



al. 1989 ). The above density can be converted to a thermal pressure in the hot medium of 

Ptherm(r) ~ lO^dynes cm" 2 (^y^j (7) 



(using a temperature of 0.67 kev; Matsumoto et al. 1997 ). This pressure could be larger than that 



required to to overcome the galaxy torque. This suggests that hydrostatic forces could dominate 
the gravitational torque. If there is a pressure gradient aligned with the jet axis in the hot medium 
then it is likely that the disk orientation could be strongly affected by a torque from the hot 
ambient medium. 



3. Modeling the Warped Disk in M84 

We first discuss constructing geometrical models for the dust features in M84 where the light 
of the galaxy is integrated taking into account the opacity from a warped disk. We then consider 
two physical models for the disk geometry which are based on a disk made up of rings of gas 
which are precessing differentially. Precession is either caused by torques from a triaxial galaxy 
or a combination of torques from an oblate galaxy and due to a pressure gradient in the ambient 
X-ray emitting medium (as proposed above). 
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3.1. A geometrical model for the warp 

We can begin with a model for the geometry of the gas disk which assumes the disk is made 



up of rings of material which precess about a given symmetry axis (e.g. used by Tubbs 1980 



Quillen ct al. 1993 , and Steiman-Cameron, Kormendy, fc Durisen 1992; ), Here we here refer to 



this axis as rh and orient it such that it is near the angular momentum axis of the gas disk. Since 
the disk rotates in a direction such that the eastern side is blueshifted with respect to the systemic 
velocity ( |Baum, Heckman, fc van Breugel 1990 ; Bower et al. 1997b ), m should be pointing near 



the north. We describe m by a position angle on the sky (3 and an inclination angle 6 which is the 
angle between this axis and the line of sight (6 = refers to a symmetry axis pointing towards us). 

We assume that a gas ring with radius r, rotation velocity v c and angular momentum 
L = rv c h precesses about rh. We define the inclination angle, u>, of a gas ring with respect to the 
symmetry axis as the angle between n and rh. The azimuthal or precession angle, a, is defined as 
the angle in the equatorial plane (plane perpendicular to rh) of the projection of n into this plane 
minus the angle of the projection of the line of sight, a increases in the direction opposite the 
rotation (e.g., see Fig. 6 of Quillen et al. 19931 ). We can then describe the geometry of the disk 



with inclination and precession angles as a function of r; 

do. 

a(r) = ao(r) = —AT + constant (8) 

u(r) = u)q. (9) 

Here we have assumed a planar initial state for the gas disk, and AT is the timescale since 
this initial condition. When the magnitude of the torque is known an estimate for the radial 
dependence of the azimuthal angle then yields an estimate for AT. When the torque is oc r _1 
then oc r _1 ; this true for the torque caused by an oblate or prolate galaxy when the rotation 
curve is nearly flat. We can then write 

a (r) = B a (^)+a c . (10) 



To produce a model we integrate the stellar light of the galaxy taking into account the opacity 
caused by the warped disk. The surface brightness profile at 1.6/iin can be well fit in the inner 20" 



with a blend of two power laws (e.g. from Ferrarese et al. 1994 ) 



I{r) = V2I C 



nY 1 

r I 



1 + 



n 



2(/?2-/3i)' 



-1/2 



(11) 



and fit as a function of semi-major axis. The break radius, r^, is where the profile changes shape 
or can also be described as near a point of maximum curvature in log log coordinates. For our 
model we assume a light density in 3 dimensions which is consistent with the above profile, 



p{r) oc s" 71 (l + s 2(72 



-71 ) 



-1/2 



(12) 
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where s = j^\Jx 2 + y 2 + z 2 /q' 2 and g is the observed isophotal axis ratio. Here we found that 



71 = 0.5 and 72 = 2 and the break radius measured at 1.6/mi (rj = 1.97 ± 0.15" « 160 pc, Quillen, 



Bower &: Stritzinger 1999] ) yields a pretty good fit to the observed surface brightness profile (with 



Pi = 0.12 db 0.04, P2 = 1-30 db 0.06). For the disk itself we assume a powerlaw form for the 
opacity of the disk (seen face on) r(r) = tq (j^m P ■ We then compared the resulting integration 
qualitatively with the morphology of the disk in an optical band (WFPC2/PC image). For the 
models shown in Fig. 2 we used tq = 1 at F547M (0.55/im) and r p = 0.2. 

Fig. 2 shows an example of such a model (referred to as Model 1) compared to the F547M 
band image. Qualitatively this model succeeds in having two self similar dust features which 
roughly correspond to the two broad features observed, one at large scales and one at smaller 
scales. The model, however does not reproduce the strong asymmetries. For example the galaxy 
dust features are wider on the eastern side than on the western side of the nucleus, but reach a 
minimum just north of the nucleus. The model on the other hand has a dust feature that is nearly 
constant width running from east of the nucleus to the north west. The observed dust features 
also extend more linearly to the west on the western side than the model does. A model which 
'wobbles' would be required to exhibit these traits. By wobbling we mean that the ring angular 
momentum vector n could follow an elliptical path (corresponding to variations in u) which has 
an azimuthally varying precession rate about this path (corresponding to azimuthal variations in 
da/dt). 



3.1.1. A Geometrical description for wobble 
We can increase the complexity of our geometric model with additional terms 

a(r) = a (r) + A a \ cos (a - ad) + A a2 sin 2(a - ad) (13) 

u){r) = lo + A,cos2(a ~ a d) (14) 

where ceo(r) comes from our more simplistic model described above, ad refers to the orientation 
angle of the maxima of the a azimuthal variations projected onto plane perpendicular to rh. 
The inclination term only contains a cos 2ao dependence because an m = 1 dependence can be 
removed with a redefinition for the symmetry axis rh. A model with these additional parameters 
(listed in Table 1, and referred to as Model 2) is also shown in Fig. 2. This model succeeds at 
reproducing the asymmetries of the observed dust features. Most of the features observed in the 
model are present in the galaxy. We infer that the dust disk is probably precessing in a way that is 
more complicated than described by the simple precession model (Eqns. 8-10) and which matches 
observations in galaxies such as Cen A ( Quillen et al. 1993j ) and NGC 4753 flSteiman-Cameron et 



al. 1992| ). We now examine what kind of physical models can produce significant extra 'wobble' 



terms in Eqns. 13 and 14. 
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3.1.2. Triaxial models 

In the above discussion we have assumed that M84 is axisymmetric, and not triaxial. To 
reproduce the alignment of dusty features we require a model with symmetry axis that is not 
aligned with the galaxy isophotal major axis. When the galaxy is triaxial, the axes of symmetry 



do not coincide with the projected isophotal axes (e.g. |de Zeeuw fc Franx 1989[) . Using equations 



for the projected angles listed in Appendix A of |de Zeeuw fe Franx (1989) and assuming a galaxy 



major axis inclination angle consistent with our fit to the warp geometry we find that either 
extreme axes ratios are required (one galaxy axis ratio lower than 0.7) or only a very small region 
(less than 10° in the azimuthal projection angle <fi of de Zeeuw Sz Franx 1989| ) of possible galaxy 



projection angles is allowed, making it unlikely. The lack of galaxy rotation (< 8km/s; Davies 



& Birkinshaw 1988) makes it unlikely that the galaxy is tumbling (as explored in |van Albada, 



Kotanyi, &: Schwarzschild 1982). 



Additional dynamical concerns also imply that the central region of M84 should not be 
triaxial. Stochasticity induced by the black hole ( |Merritt &: Valuri 1996 ) should cause triaxiality 



to be short lived (of order the local dynamical time) in the central regions of the galaxy. There 



is little twist in the isophotes of M84 ( van den Bosch et al. 1994 , Peletier et al. 1990 ) which 



should be observed if the galaxy is strongly triaxial and not oriented with an axis coincident with 
the line of sight. We also note that the warped disk appearance and its misalignment with the 
galaxy major axis persists to well within the break radius ~ 2" ~ 160 pc of the nucleus where the 
ellipticity is slightly reduced and the slight boxiness disappears (as seen in the NICMOS 1.6/um 
isophotes; |Quillcn, Bower fc Stritzinger 1999 ). Mechanisms that reduce ellipticity and boxiness 



(such as scattering) are likely to destroy triaxiality as well. In addition van Dokkum & Franx 



1995) found that triaxiality could not explain the observed disk/galaxy misalignments in a sample 
of elliptical galaxies. Whereas precession times are similar, disk settling times in triaxial galaxies 
are substantially faster than in axisymmetric systems ( |Habe &: Ikeuchi 1985 ). 



We can also determine if a triaxial model is viable based on the size and form of the 'wobble' 
terms introduced above, which are probably needed to describe the disk morphology. We associate 
the symmetry axis of the disk m with one symmetry axis of the triaxial galaxy. The amplitude 
of variation in inclination angle described by A w can be estimated from the azimuthal variation 
in the torque on a gas ring as the ring precesses. This results in a torque that causes the ring to 
change inclination (lo) rather than precession rate. This amplitude should be roughly equivalent 
to the size of the ellipticity of the gravitational potential in the plane perpendicular to fh times the 
angle to (for to small). Since this angle is small in M84 and the ellipticity of the potential should 
be ~ 1/3 the axis ratio of the density in this plane, should be quite small (we estimate less 
than a degree). This is far less than exhibited in our Model 2 described above. In a triaxial galaxy 
it makes sense that this amplitude should be small because the gravitational potential is always 
smoother than the actual density distribution. For a triaxial galaxy we also expect no sin(ao) 
variation in da/dt (or A a i term) which were a part of Model 2 discussed above. Dynamical models 
which neglect triaxiality for the disk kinematics and morphology in Centaurus A (e.g. {Bparkcj 
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1996, Quillen et al. 1993 ) are successful despite the evidence that this galaxy probably is triaxial 
(Hui et al. 1995). This supports our statement that triaxial galaxies are unlikely to result in disk 
precession models with large 'wobble' terms. 



3.2. Warp model including a pressure torque 

Above we proposed that the ambient X-ray emitting ISM could exert sufficient torque on the 
gas disk to affect the orientation of the gas disk. Here we elaborate on this possibility and search 
for a model that can accurately match the inferred geometry of the gas disk. 

The disk rotates in a direction such that the eastern side is blueshifted with respect to the 
systemic velocity (Baum et al. (1990); Bower et al. 1997b| ) . The gas disk, if warped, has orientation 



such that with increasing radius the azimuthal angle of greatest disk inclination moves counter 
(retrograde) to the direction of rotation. This orientation would be consistent with an oblate 



galaxy potential (if the warp were indeed caused by the torque from the galaxy; e.g., |Tubbs 1980[) . 
For an oblate galaxy potential the shape of the potential causes a force towards the galactic plane 
of symmetry, resulting in a 'retrograde' warp. If higher pressures exist in the ISM along the jet 
axis then the pressure gradient in this gas would exert a force on the gas disk towards the plane 
perpendicular to the jets. When the precession period decreases with increasing radius this again 
would cause a 'retrograde' warp but about a symmetry axis aligned with the jet rather than a 
galactic axis of symmetry. 

We can estimate the torque resulting from the hot ISM as follows: As typically assumed for a 
non-axisymmetric galaxy, we can describe the pressure in the volume filling X-ray emitting ISM as 
having an the axis of symmetry, p, of the isobars in the hot gas. We chose to orient p with a sign 
such that it is nearest the average angular momentum vector of the disk. The pressure gradient 
across a gas disk with mass surface density £ and thickness h (corresponding to the vertical scale 
height) results in an average torque per unit mass (around a ring of radius r) 

h dP 

^ = 2s m (15) 

where the direction of the torque on the ring is given by p x h. 

The average torque per unit mass from a non-axisymmetric gravitational potential can be 
determined similarly with 

f g = e$>v 2 c cos 9 g {n x g) (16) 

where we have assumed a cos(2(9 9 ) form for the non-axisymmetric part of the gravitational 
potential. Here g refers to the ring inclination with respect to the axis of symmetry of the galaxy, 
or the angle between n and g, for g the galaxy symmetry axis and e$ is the ellipticity of the 
gravitational potential. The total average torque per unit mass on a gas ring would then be given 
by the sum of the two torques f = fg + f p . 
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Precession of the gas disk takes place approximately about a vector 

- dP h - 2~ 

m = W2S P + W (17) 

where Q£ is estimated at the average angle between h and p. We can define orientation angles 
for the gas disk as a function of radius with respect to this vector rh. As in §3.1 we define the 
inclination angle, uj, and a for a ring with angular momentum axis n with respect to the precession 
axis rh. The total torque is 

r rj rh x n (18) 

and the angular precession rate about rh is 

da rh - h 



dt rv r 



(19) 



When the galaxy is oblate the gas disk will precess about an axis intermediate between the 
galaxy axis of symmetry {g) and the pressure axis of symmetry (p, which we associate with the jet 
axis). However when the galaxy is prolate the sign of the galactic torque term changes (equivalent 
to setting e$ to be negative) and precession could occur about an axis which is not obviously 
between the two symmetry axes. 

Because of the sense of the warp, in M84 we expect that the galaxy is nearly oblate and so 
that the gas disk should precess about an axis intermediate between the two symmetry axes. In 
fact we can measure the orientation of this axis and from it estimate the size of the pressure torque 
compared to the the galactic term (assuming that the galaxy is not strongly triaxial). 



3.2.1. Ratio of pressure to galactic torque 

In M84 we observe that the dust features are aligned at an angle (PA= 65 — > 80°) intermediate 
between the jet axis (PA= 0; Jones et al. 1981 ) and the galaxy major axis (PA= 129°, |van Dokkuml 



fc Franx 1995| ). Because the north and south jets are have similar intensities we assume that the 



jet axis is perpendicular to the line of sight. We find a rough fit to the warp with a symmetry axis 
with PA between —10 and —15°. So the angle between the jet and rh is about 10 — 15° and the 
angle between rh and the galaxy axis is about 35 — 40°. This suggests that the torque from the 
pressure is larger than the galactic torque so that 

^ « 2 - 4 (20) 

Tp 

This lets us estimate the size of the azimuthal pressure component. For ^ = e p Po(r) 
Po(r) ~ 3^ (21) 



e p h 



3 x 10 dynes cm 



-1 / „, \2 / v \ / ^ \ -1 



e$ \ / e p \ / v c \ / E \ fh/r 



0.05/ V 1.00/ V300kms _i / V2M p C - 2 / V 0.1 / UOOpc 
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Here e p represents the ellipticity of the isobars. We have adopted a range for the surface density £ 
based on the extinction in the dust features, van Dokkum & Franx (1995) and Bower et al. (1997a) 



estimate a total disk mass of 1O 6 M0 and 9 x 10 6 Mq respectively. An estimate for the disk mass 
based on the NICMOS images is consistent with the low end ~ 1O 6 M0. For a disk of constant 
surface density (consistent with a small r p , see §3.1) truncated at r = 400 pc the lower of these 
mass estimates give £ ~ 2M pc~ 2 . 

Although there are uncertainties in the parameters, the pressure we estimate above is 
similar to the thermal pressure we estimated from the ambient volume filling X-ray emission (see 
Eqn. 7). The alignment and precession of the disk therefore could be consistent with the model 
proposed here where a pressure exerted by the X-ray medium (aligned with the jet) and the 
non-axisymmetric galaxy both exert torques on the gas disk. 



3.2.2. Warp timescale 

As mentioned in §3.1, once the size and radial form of the torque is known, the radial 
dependence of a can be used to estimate a timescale. Here we have assumed a planar initial state 
for the gas disk, and AT is the timescale since this initial condition (see Eqn. 2). For a warp 
caused purely by torque from the galaxy da/dt ~ e$f2 oc r~ l when the rotation curve is nearly 
flat. From our modeling we estimate B a = 1.2 ± 0.3 x 10 3 radians pc yielding an estimate of 

AT = 1.8 x 10 V ( ^ ) (-^-) ^ ( VC r) 1 1 (22) 

y \1.2 x 10 3 radians pcAo.05/ \ 300km s -1 / V 3 J K J 

This suggests that precession has occurred only about 10 7 years since an event responsible 
for distributing the gas and dust, or causing the present orientation of the jets. Our model places 
time limits on the stability of the jet orientation. Here we link the morphology of the disk to the 
history of pressure asymmetries caused by the jets themselves. If the jets were at a different angle 
previously (as suggested from the large scale radio emission) then the large scale morphology of 
the disk could be linked to the past jet orientation. An alternate mechanism is then required to 
change the jet orientation in a way not related to the disk orientation (for example as proposed 
by [Roos, Kaastra, fc Hummel 1995 involving binary black holes). 



3.2.3. Wobble 

In general the sum of the two torques will lead to an elliptical path (or 'wobble') for the 
angular momentum vector of a gas ring, n, about an axis ~ rh which also varies in angular 
precession rate. This wobble corresponds to the variation in torque that a ring feels as it precesses. 
From a measurement of this 'wobble' we can determine if the pressure gradient is strongly 
dependent on 6 p , the angle between h and p. As in the case of a triaxial galaxy, for a cos 28 p 
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dependence we would expect that the wobble about simple precession should be small (less than 
a degree). However we achieve a better correspondence to the observed morphology with larger 
azimuthal variations in uj and da/dt suggesting that the pressure is strongly dependent on the 
inclination angle p from the jet axis. 

We now estimate the difference in the torque as the angular momentum axis precessed about 
fh As the ring precesses variation in the torque reach extremes nearest and furthest from the jet 
axis of the size 



This results in a variation in that is proportional to cos(«o) (described in our model by A a i). 
The component of At in the plane perpendicular to fh causes the ring to change inclination and 
is about the same size as At times an additional factor of uj and results in a cos(2ao) dependence 
of uj which we described in our model by A^. (The cos(ao) dependence can be removed by a 
redefinition of fh). Here our estimates have assumed that uj is small. 

If both the pressure and gravitational torques show the same angular dependence (e.g. 
cos(2#)) then the dependence of a cancels to first order in the angles 9 g and 9 P and ^ varies only 
to second order in uj. The dependence on a can be approximated as 



where m refers to the angular dependence of ^ on 9 P and the variation in uj is ~ A a \uj. a(r) 
is then described by Eqns. 8, 10 and 13 and uj(r) is described by Eqn. 14. The precession rate is 
expected to be fastest when the ring angular momentum vector is furthest from the pressure or 
jet axis, p, and nearest g, the galaxy symmetry axis. From this we can see that large azimuthal 
variations in the angular precession rate and inclination could be symptoms of a strong angular 
dependence of the isobars. The signs of A w and A a \ in our Model 2 (see Table 1) suggests that 
the galaxy axis (to the north-west on the sky) is pointing somewhat towards us with the angle 
between it (projected on to the planet perpendicular to fh) and the line of sight |a^| ~ 50°. 



We have ignored the possibility that the isobars might actually be affected by the galaxy 
potential. This situation would arise naturally to some extent because as energy is lost from 
the jet into the ISM energy would be transfered in the direction of least gravitational force, i.e. 
towards the galactic major axis. However if hydrostatic equilibrium applies, the sound speed of 
the gas sets a particular scale length in the hot gas whereby 7 ~ f^- Because the sound speed 
in the hot gas c s ~ v c ~ 300 km/s hydrostatic equilibrium would set h/r ~ 1 and so the isobars 
are unlikely to be significantly affected by the gravitational potential, and are more likely to be 
determined by the nature of the dissipation from the jet into the ISM. 




(23) 



A a \ m uJo(m 2 — 4) cos(«o) 



(24) 



3.24. A note 
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4. Summary and Discussion 

In this paper we have estimated the timescale of a gas disk to precess in the non-spherical 
gravitational potential of M84. This timescale is a few times 10 years at 100 pc where the 
dust features are misaligned with the galaxy isophotes. For the disk to remain misaligned with 
the galaxy potential some mechanism must operate faster than this. While a cooling flow could 
replenish the disk on this timescale it is difficult to explain why the disk is at a roughly constant 
angle within r < 7" and yet twists at this radius forming an overall S-shape in the Ha + [Nil] 
emission ( |Baum et al. 1988 ). Extremely fast accretion through the disk itself would require 



a substantial gas reservoir at large radii or a very short disk lifetime. It would also require a 
place to put the excess accreted gas mass, such as a wind, an inner disk or advection-dominated 
accretion into the black hole. A combination of fast accretion and replenishment by a cooling flow 
could possibly result in inflow faster than the precession rate, but it is not clear whether this 
combination could account for the disk morphology. The AGN is not luminous enough for the 



radiative induced warp mechanism of Pringle (1996) to operate. None of these possibilities would 



provide a good explanation for the observed morphology of the gas and dust disk. 

As an alternative to these external mechanisms we consider the possibility of a local force on 
the disk. We estimate the pressure required to overcome the torque from the galaxy and find that 
it is small compared to the thermal pressure inferred from X-ray observations. Pressure gradients 
in this ambient hot ISM could therefore overcome the galaxy torque. We therefore propose that 
pressure gradients in an energetic low density medium in M84 could strongly affect the orientation 
of the gas disk on the scale of a few hundred pc. 

By integrating the light of the galaxy through a dusty warped disk we find that the gas disk 
in M84 is likely to differ from a simple precession model where the precession rate is constant 
with azimuthal angle and the angular momentum axis of a gas ring traces a circular path. A 
triaxial model for the galaxy, though it would explain the misalignment of the dust features with 
the galaxy isophotal major axis, is not a good explanation for a variety of reasons. Because of 
the expected weak dependence of the gravitational torque on the position angle wobble of the 
angular momentum axis of a gas ring is likely to be smaller than that we infer from the disk 
geometry. The misalignment also persists to small radii where the galaxy is expected to be oblate 
and not triaxial. We propose instead that the morphology of the gas disk in M84 is consistent 
with a warped geometry where precession is caused by a combination of a galactic torque and a 
torque due to pressure gradients in the ambient X-ray emitting gas. The alignment of the disk can 
be used to estimate the ratio of these torques. Precession occurs at an axis between the jet and 
galaxy major axis, but nearer to the jet axis. This angle implies that the pressure torque is 2-4 
times larger than the galactic torque. A rough model to the morphology of gas disk also allows 
us to estimate the degree of precession that has taken place. Assuming the initial condition of 
gas in a plane we estimate the timescale since then to be a few times 10 years. A better model 
to the morphology of the disk is achieved when precession takes place about an elliptical rather 
than circular path and precesses fastest when the angular momentum vector is furthest from the 
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jet axis. This suggests that the isobars are strongly dependent on angle from the jet axis. 

Recent investigations find that almost all ellipticals have dust ( |van den Bosch et al. 1994 ). 
In cases of non-active elliptical galaxies having dusty disks misaligned with the galaxy isophotes, 
the misalignment must be caused by another mechanism (perhaps a cooling flow or a galaxy 
merger), van Dokkum & Frame's (1995) sample contains only five galaxies with radio power 
logP„ (6 cm)< 20 (W/Hz) and misaligned dust features. Further investigation is needed with a 
larger sample to determine how common these non-AGN misaligned cases are and if the models 
proposed here are applicable. One consequence of our proposed mechanism is that on very long 
timescales we expect the disk to become multiply warped or rippled. On a timescale of roughly 10 
times the precession time we would then expect the disk to settle into a quasi-stationary surface 
nearly perpendicular to the jet axis. This might account for the observed alignments between jets 
and dust features. Comparison of high resolution X-ray morphology (such as will be possible with 
AXAF), dust morphologies and ionized gas kinematics should determine if the type of models 
introduced here are appropriate. 

We acknowledge helpful discussions and correspondence with E. Emsellem, A. Eckart, 
R. Green, G. Rieke, M. Rieke, G. Schmidt, D. Hines, P. Pinto, D. DeYoung and F. Melia. 
Support for this work was provided by NASA through grant number GO-07868.01-96A from the 
Space Telescope Institute, which is operated by the Association of Universities for Research in 
Astronomy, Incorporated, under NASA contract NAS5-26555. ACQ also acknowledges support 
from NSF grant AST-9529190 to M. and G. Rieke and NASA project no. NAG-53359. GB 
acknowledges support from the STIS Investigation Definition Team. 



-15- 



REFERENCES 

Abramowicz, M. A., Czerny, B., Lasota, J.-P., & Szuszkiewicz, E. 1988, ApJ, 332, 646 

Abramowicz, M. A., & Lasota, J.-P. 1996, Comments Astrophysics, 18, 141 

Axon, D. J., Unger, S. W., Pedlar, A., Meurs, E. J. A., Whittle, D. M. & Ward, M. J. 1989, 
Nature, 341, 631 

Bardeen, J. M., & Petterson, J. A. 1975, ApJ, 193, L65 

Baum, S. A. Heckman, T., Bridle, A, van Breugel, W., & Miley, G. 1988, ApJS, 68, 643 
Baum, S. A., Heckman, T., & van Breugel, W. 1990, ApJS, 74, 389 
Baum, S. A., Heckman, T. M., & van Breugel, W. 1992, ApJ, 389, 208 

Begelman, M. C. 1982, in Astrophysical Jets: Proceedings of an international workshphelp in 

Torino, Italy 1982, edited by A. Ferrari and A. G. Pacholczyk, (Reidel, Dordrecht, Holland), 
p 215. 

Binney, J., & Tremaine, S. 1987, Galactic Dynamics, (Princeton U. Press). 

Bower, G. A., Heckman, T. M., Wilson, A. S., & Richstone, D. O. 1997a, ApJ, 483, L33 

Bower, G. A., et al. 1997b, ApJ, 472, Llll 

Cecil, G., Morse, J. A., & Veilleux, S. 1995, ApJ, 452, 613 

Davies, R. L., & Birkinshaw, M. 1988, ApJS, 68, 409 

de Zeeuw, T., & Franx, M. 1989, ApJ, 343, 617 

De Young, D. S. 1981, Nature, 293, 43 

Ferrarese, L., van den Bosch, F. C, Ford, H. C, Jaffe, W., & O'Connell, R. W. 1994, AJ, 108, 
1598 

Ferrarese, L., Ford, H. C, & Jaffe W. 1996, ApJ, 470, 444 

Gebhardt, K., Richstone, D., Ajhar, E. A., Lauer, T. R., Byun, Y.-L, Kormendy, J., Dressier, A., 
Faber, S. M., Grillmair, C, & Tremaine, S. 1996, AJ, 112, 105 

Gunn, J. E. 1979, in Active Galactic Nuclei, edited by C. Hazard and S. Mitton, (Cambridge 
University Press, Cambridge), p. 213 

Habe, A., & Ikeuchi, S. 1985, ApJ, 289, 540 

Hansen, L., Norgaard-Nielsen, H. U., & Jorgensen, H. E. 1985, A&A, 149, 442 



-16- 



Harms, R. J., et al. 1994, ApJ, 435, L35 

Heckman, T. M., Baum, S. A., van Breugel, W. J. M, & McCarthy, P. 1989, ApJ, 338, 48 

Herbig, T., & Readhead, A. C. S. 1992, ApJS, 81, 83 

Hernquist, L., Heyl, J. S., & Spergel, D. N. 1993, ApJ, 416, L9 

Hui, X., Ford, H. C, Freeman, K. C, & Dopita, M. A. 1995, ApJ, 449, 592 

Jaffe, W., Ford, H. C, O'Connell, R. W., & van den Bosch, F. C, & Ferrarese, L. 1994, AJ, 108, 
1567 

Jones, D. L., Sramek, R. A., & Terzian, Y. 1981, ApJ, 246, 28 
Kotanyi, C, & Ekers, R. 1979, A&A, 73, LI 

Maloney, P. R., Begelman, M. C, & Pringle, J. E. 1996, ApJ, 472, 582 

Matsumoto, H., Koyama, K., Awaki, H., Tsuru, T., Loewenstein, M., & Matsushita, K. 1997, 
ApJ, 482, 133 

Merritt, D., & Valuri, M 1996, ApJ, 471, 82 

Mollenhoff, C, Hummel, E., & Bender, R. 1992, A&A, 255, 35 

Mould, J., et al. 1995, ApJ, 449, 413 

Pelletier, R. F., Davies, R. L., Illingworth, G. D., Davis, L. E., & Cawson, M. 1990, AJ, 100, 1091 
Pringle, J. E. 1996, MNRAS, 281, 357 

Quillen, A. C, Bower, G. A., & Stritzinger, M. in preparation 

Quillen, A. C, de Zeeuw, P. T., Phinney, E. S., & Phillips, T. G. 1992, ApJ, 391, 121 
Quillen, A. C, Graham, J. R., & Frogel, J. A. 1993, ApJ, 412, 550 
Rees, M. J. 1978, Nature, 275, 516 
Rees, M. J. 1984, ARA&A, 22, 271 

Rees, M. J., Begelman, M. C, Blandford, R. D., & Phinney, E. S. 1982, Nature, 295, 17 
Roos, N., Kaastra, J. S, & Hummel, C. A. 1993, ApJ, 409, 130 
Sansom, A. E., et al. 1987, MNRAS, 229, 15 
Sparke, L. S. 1996, ApJ, 473, 810 



-17- 



Steiman-Cameron, T. Y., & Durisen, R. H. 1988, ApJ, 325, 26 

Steiman-Cameron, T. Y., Kormendy, J. & Durisen, R. H. 1992, AJ, 104, 1339 

Tadhunter, C. N., Fosbury, R. A. E., & Quinn, P. J. 1989, MNRAS, 240, 225 

Thomas, P. A., Fabian, A. C, Arnaud, K. A., Forman, W., & Jones, C. 1986, MNRAS, 222, 655 

Tubbs, A. D. 1980, ApJ, 241, 969 

van Albada, T. S., Kotanyi, C. G., & Schwarzschild, M. 1982, MNRAS, 198, 303 

van den Bosch, F. C, Ferrarese, L., Jaffe, W., Ford, H. C, & O'Connell, R. W. 1994, AJ, 108, 
1579 

van Dokkum, P. D., & Franx, M. 1995, AJ, 110, 2027 



-18- 



This preprint was prepared with the AAS IAT^X macros v4.0. 



-19- 



Fig. 1. — The relative orientation of the gas disk and stellar isophotes in M84. HST/WFPC2/PC 
continuum imaging of the center of M84 from Bower et al. (1997a), with a resolution of 0.1" (8 
pc). The grayscale shows the (V-I) map, with values of (V-I) ranging linearly from 2.0 (black) to 
1.3 (white). The reddest color in the data is (V-I) = 1.7. Contours from the V-band image are 
superimposed, where the contour interval corresponds to a factor of \[2 in intensity. The jet has 
PA f» 10° on the pc scale and PA ps 0° on the 100 pc scale (Jones, Sramek, Sz Terzian 1981). The 
dotted lines show the approximate morphology of the S-shape twist in the larger scale Ha + [Nil] 
emission map of Baum et al. (1988). Within these dotted lines the large scale emission is alligned 
with the dust lanes. North is up and East is to the left. 



Fig. 2. — Models compared to the WFPC2/PC F547M image. These figures are oriented with 
precession symmetry axis of the disk oriented up. North is at a position angle 10° east of the y 
axis, a) A circular precession model with no wobble (Model 1 with parameters listed in Table 1). 
b) A model with wobble (Model 2). c) F547M image. 
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Table 1. Model Parameters for M84's Disk 



Model 


P 


e 


B a 


a c 


LOO 


A al 




Ctd 




deg 


deg 


radians pc 


deg 


deg 


deg 


deg 


deg 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


-10 


112 


-1.2 x 10 3 


180 


12 











2 


-10 


112 


-1.2 x 10 3 


180 


8 


-20 


8 


-50 



NOTES.- Parameters are discussed in §3.1. Columns: (1) Model 1 has no 'wobble' and is shown in Fig. 2a. Model 
2 is shown in Fig. 2b; (2) Position angle on the sky (E of N) of symmetry axis of warp, m; (3) Inclination of m. 
(6 = refers to a symmetry axis pointing towards us); (4) Dependence of ao on radius, See Eqns. 8, 10 and 13; (5) 
Constant term for qq. See Eqns. 10 and 13; (6) Average disk inclination with respect to the warp symmetry axis, 
(angle between rh and n); see Eqns. 9 and 14. (7) Amplitude of m = 1 azimuthal variation in a, see Eqn. 13; (8) 
Amplitude of m = 2 azimuthal variation in u, see Eqn. 14; (9) Angle where extrema of a and lo variation occurs. See 
Eqns. 13 and 14. 
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